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ABSTRACT 



Context. Polarimetry of extrasolar planets becomes a new tool for their investigation, which requires the development of diagnostic 
techniques and parameter case studies. 

Aims. Our goal is to develop a theoretical model which can be applied to interpret polarimetric observations of extrasolar planets. 
Here we present a theoretical parameter study that shows the influence of the various involved parameters on the polarization curves. 
Furthermore, we investigate the robustness of the fitting procedure. We focus on the diagnostics of orbital parameters and the estima- 
tion of the scattering radius of the planet. 

Methods. We employ the physics of Rayleigh scattering to obtain polarization curves of an unresolved extrasolar planet. Calculations 
are made for two cases: (i) assuming an angular distribution for the intensity of the scattered light as from a Lambert sphere and 
for polarization as from a Rayleigh-type scatterer, and (ii) assuming that both the intensity and polarization of the scattered light are 
distributed according to the Rayleigh law. We show that the difference between these two cases is negligible for the shapes of the 
polarization curves. In addition, we take the size of the host star into account, which is relevant for hot Jupiters orbiting giant stars. 
Results. We discuss the influence of the inclination of the planetary orbit, the position angle of the ascending node, and the eccen- 
tricity on the linearly polarized light curves both in Stokes Q/I and U/I. We also analyze errors that arise from the assumption of a 
point-like star in numerical modeling of polarization as compared to consistent calculations accounting for the finite size of the host 
star We find that errors due to the point-like star approximation are reduced with the size of the orbit, but still amount to about 5% for 
known hot lupiters. Recovering orbital parameters from simulated data is shown to be very robust even for very noisy data because 
the polarization curves react sensitively to changes in the shape and orientation of the orbit. 

Conclusions. The proposed model successfully diagnoses orbital parameters of extrasolar planets and can also be applied to predict 
polarization curves of known exoplanets. Polarization curves of extrasolar planets thus provide an ideal tool to determine parameters 
that are difficult to obtain with other methods, namely inclination and position angle of the ascending node of orbits as well as true 
masses of extrasolar planets. 

Key words, polarization ~ scattering - methods: numerical - eclipses - stars: planetary systems 



1. Introduction 

Direct observations of extrasolar planets succeeded only re- 
cently by measuring infrared thermal emission with the Spitzer 
, Space Telescope (Charbonneau et al. 2005; Deming et al. 2005, 
i l2006t [H arrington et al. 2006), polarimetric measurements 
' (iBerdvu'gina et al., i2008.) . and direct imaging (iMarois et al.l 
' 120081) . Polarimetry as a method to investigate extrasolar plan- 
ets is particularly valuable because it involves direct observation 
of light scattered by the planetary atmosphere and reveals infor- 
mation about the chemical composition, the physical structure, 
and the evolution of the planetary atmosphere. So far the iden- 
tification of atmospheric constituents is limi ted to transmission 
spectroscopy of transiting extrasolar planets ("Charbonne au et al] 
[2002; Vidal-Madiaret al..,2003, 2004; Tinetti et al. 2007) . while 
polarimetry is capable of expanding these studies for non- 
transiting planets. 

The fi rst successful pol a rimetr ic detection of an extrasolar 
planet bv iBerdvugina et al] ( 20081) (see also new da ta and dis- 
cussion on non-detections in Berdvugina et alJl2009l) opens new 
perspectives for diagnosing the physical properties of planetary 
atmospheres and underlines the urgent need for the development 
of theoretical tools to interpret the data. Polarization signals of 



extrasolar planets have akeady been numerically simulated by 
several authors, who in particular investigated the influence of 
various scatter ing particles r esulting from different atmosp heric 
compositions. Seager e t"al] (l2000l) . ISaar & Seageii ( 12003 ), and 
Hough & Lucas ( 2003) have discussed close-in extrasolar giant 
planets with orbital radii smaller than 0.05 AU, which are thus 
unresolved and relatively hot, while iStam et al.l (|2004|) studied 
the flux and polarized spectra of resolved Jupiter-like pla nets in 
larger orbits. In contrast to others Sengupta & Maitil (12006) have 
considered elliptical rather than circular orbits to calculate the 
polarization in the R passband due to scattering by water and sil- 
icate condensates, taking the planetary oblateness also into ac- 
count. 

Most of these previous numerical simulations present light 
curves or spectra of the total degree of polarization only. They 
additionally are restricted by the approximation of the host stars 
as point-sources. However, the actual observables obtained from 
polarimeters are the Stokes parameters, not the total degree of 
polarization. It is obvious that part of the information concern- 
ing the orbital parameters is lost when forming the total degree 
of polarization. Moreover, the polarimetric method is well suited 
to detect extrasolar planets that closely orbit giant stars where the 
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description of the host star as a point-source becomes question- 
able. 

The numerical work has mainly concentrated so far on the 
red and near infrared parts of the spectrum, where the detailed 
composition of the atmosphere and the shape and size of the 
scattering particles play an important role. Also, the unsuccess- 
ful polarimetric detectio ns of ex oplanets were reported for red 
and green wavelengths lLucas_et al. 2009; Wiktorowicz 2009). 
The first successful detection of polarization from an extra- 
solar planet was obtained in the B and U passbands though 
(iBerdvupna et al. l2008. 2009) f or HD18 9733b. As thoroughly 
discussed and modeled by Berdvugina et al. (2009), this clearly 
indicates an atmosphere in which Rayleigh scattering plays a 
dominant role. The same conclusion was re ached in a multi 
wav elength analysis of planetary transit curves (iPont et al.l2008t 
iLec avelier Des Etangs et al. 2008). In this atmosphere, a max- 
imum polarization amplitude in green (red) wavelengths is at 
least a factor of two (ten) smaller than that in the B-band. This 
would explain the non-detections by iLucas et aT 1 (l2009l) and 
IWiktorowicj (1200 9^. In this paper we concentrate on Rayleigh 
scattering and exclude Mie scattering, as the latter has not yet 
been detected in any extrasolar planet. We leave the case of Mie 
scattering for later studies. 

We discuss numerical simulations of polarized light curves 
obtained from Rayleigh scattering on unresolved extrasolar plan- 
ets, particularly actually observed Stokes Qfl and Ufl phase 
curves. The modeling strategy introduced here has been ap- 
plied for the interpretation of the recently obtained polarimetric 
data for HD 189733b (Berdyugina et al. 2008). We describe here 
the numerical model in detail, perform a theoretical parameter 
study for involved parameters, verify the validity of the Lambert 
sphere approximation, and test the stability of the procedure to 
fit observations. At this stage we concentrate on orbital parame- 
ters, in general for eccentric orbits, and on the planetary radius, 
which represents the first step when interpreting data and is suf- 
ficient for the currently available observations and polarimetric 
precision. Particular attention is paid to the correct treatment of 
the size of the host star. 



2. Modeling polarization of extrasolar planets 

2. 1 . Model assumptions 

In this section we describe the numerical model for computing 
the polarization signal from extrasolar planets. We employ the 
following main assumptions: 

- the intensity of the scattered light is given either by the 
Lambert sphere or by Rayleigh scattering, 

- the polarization is defined by the Rayleigh scattering law, 

- the host star has a finite size, 

- the direct irradiance from the star is unpolarized, 

- the planet is unresolved from the host star, 

- the radius of the planet is much smaller than the stellar ra- 
dius, 

- the thermal emission from the planet is neglected. 

The first two assumptions form two cases: (i) a combination 
of the Lambert sphere angular distribution for intensity of the 
scattered light and the Rayleigh phase function for polarization 
and (ii) scattering consistently described by the Rayleigh law for 
both intensity and polarization. The first case was employed pre- 
viously by several authors as it has an analytical solution. Here 
we compare this solution with the Rayleigh scattering calcula- 



tion and demonstrate that both cases result in very similar shapes 
of the polarization curves. 

By considering the finite size of the host star we have to take 
into account that different parts on the stellar surface have in 
general individual distances to the planet because of the spher- 
ical shape of the star. In our calculations the surface of the star 
is divided into individual pixels. They are chosen to be small 
enough so that all light rays emitted from one stellar surface el- 
ement arrive nearly parallel on the planetary atmosphere, which 
requires our assumption that the planet is small compared to the 
star. Therefore the incident radiation from one single stellar pixel 
results in a contribution to the scattered radiation that can be 
computed as in the case of a point-like star In general, the in- 
cident flux, the scattering angle, and the phase angle differ for 
every pixel. The total scattered radiation field is obtained simply 
by adding up the contributions from all stellar surface elements. 
It is crucial to consider only the actually visible stellar surface as 
seen from the position of the planet, which approaches exactly 
one stellar hemisphere for very distant planets, but is in general 
smaller than that for close-by planets. 

The ideal size of the stellar pixel in principle depends on 
the stellar radius and the size of the planetary orbit. For the cal- 
culations presented in this paper we chose a pixel size of one 
degree both in longitude and latitude. Pixel sizes of more than 
five degrees have occasionally lead to inaccurate results in test 
calculations. For smaller pixels the computed flux and polariza- 
tion were always independent of the pixel size, as long as it is 
still high enough to prevent numerical errors due to the machine 
precision. At the end of this section we briefly review the re- 
sults for a point-like host star compared with an extended star 
for completeness' sake. 

2.2. Lambert-Rayleigh approximation 

The flux F from the unresolved planetary system measured by a 
distant observer consists of two parts, namely 

F = F^+F,,, (1) 

where is the direct stellar flux and F^c represents the frac- 
tion of the stellar flux which is scattered by the extrasolar planet 
into the line-of-sight. We neglected the thermal emission from 
the planet because we are mainly interested in the optical wave- 
lengths at this stage. 

The direct flux is given by 

F^ = 2n-^ Ifidfi , (2) 

with R-t. the stellar radius, d the distance from the star to the 
observer that is assumed to be very large compared to the semi- 
major axis of the planetary orbit, / the intensity at the stellar 
surface, and ju = cos 6, 6 being the angle between the line-of- 
sight and a stellar surface element as seen from the center of the 
star. The integral accounts for limb darkening too. The absolute 
value of the flux is of no concern for the polarization, which is 
a relative quantity. We define it by choosing the normalization 
condition 

IttRI f Ifidn = 1 , (3) 
Jo 
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which scales the intensity at the surface of the star. The direct 
stellar flux then simplifies to = We have modeled the 
limb darkening with the expression introduced bv iClaret (,2000i) 



(4) 



k=i 



In this paper we used two sets of coefficients, one for a hypothet- 
ical exoplanet orbiting a giant star (oi - 0.4908, a2 = -0.5346, 
aj, = 1.2703, fl4 = -0.2889) and another one for the Sun (oi - 
0.4767, a2 = -0. 1 59 1 , fl3 = 1 -07 1 1 , fl4 = -0.5 1 54), bo th approxi- 
mating the limb darkening in the B passband ( IClaret ''2000). 

To obtain the scattered flux Fsc we considered the total en- 
ergy incident on the planet per second and frequency unit from 
one single stellar surface element 



(5) 



where the index "p" indicates that the corresponding quantity 
depends on the selected stellar surface element or pixel. The pa- 
rameter r represents the radius of the planet, i.e. nr^ is just the 
cross-section of the planet, /p is the intensity emitted from the 
surface element in the direction of the planet, and dQ^ is the solid 
angle covered by the stellar surface element as seen from the 
planet and should not be confused with the position angle of the 
ascending node Q.. The distance from the stellar surface element 
to the planet is implicitly contained within the solid angle di3p 
and thus dilTers for different pixels e.g. at the limb or disk center. 
If the planet is small compared to the star, we can assume that 
d£?p coincides for all positions on the planet. Correspondingly, 
we compute AQ^ always for the center of the mass of the planet. 
From the incident light Fi,ic,p a fraction 

/scp = -P c[>(ap) (6) 
n 

is scattered towards the observer. Here I/tt represents the frac- 
tion of light scattered back to the star by a Lambert disk at phase 
angle zero, p is the geometric albedo, which is 2/3 for a Lambert 
sphere, and a-p is the phase angle, i.e. the angle between a stellar 
surface element and the observer as seen from the planet. For a 
Lambert spher e the phase function (J takes the analytical form 
(iRussellll 19161) 



<I'(ap) = — (sin a-p + {n — o-p) cos o-p) 



(7) 



Combining Eqs. (|5]i and ^ we find the contribution to the flux 
that is scattered by the planet towards a distant observer 



" «7 visible 



(8) 



The integral is carried out over the visible surface of the star as 
seen from the planet. It is easily verified that the scattered flux 
Fsc as given in Eq. ([8]) approaches the simplified case of the star 
as a point-source for large planetary orbits, i.e. the second part 
given in Eq. ( fT6] l. 

Scattering polarization was added in an ad hoc fashion, be- 
cause a Lambert sphere scatters the incident radiation isotrop- 
ically and in principle unpolarized. We estimated t he po lariza- 
tion in a similar way as was done by ISeaeer et al.l (I200Q) for a 
Lambert sphere, thus creating an idealized case, which is sufli- 
cient to discuss the influence of orbital parameters on the polar- 
ization curves. The calculated degree of polarization defines an 
upper limit, so that an inferred planetary radius corresponds to a 
lower limit. 



We employed two different reference frames to represent the 
Stokes Q and U parameters. The final results are all given with 
respect to celestial north for the observer on Earth, so that a pos- 
itive Stokes Q describes light that is linearly polarized parallel 
to the local meridian of the observer. This first reference frame 
is referred to as the "observer's frame". A different reference 
frame is temporarily used while adding up contributions from 
different stellar surface elements, where we define the Stokes 
parameters with respect to the scattering plane so that positive 
Stokes Q represents light that is linearly polarized perpendicular 
to the scattering plane. 

Here we assumed that light incident on the planet is scattered 
and that the resulting Stokes Q and U are defined by Rayleigh 
scattering. Considering scattered light originating from one stel- 
lar surface element, and choosing the scattering plane as the ref- 
erence frame (indicated by primed variables), we obtained a con- 
tribution to Stokes Q' given by 



1 J pr J 

6p = ^^'incp/scpsm lAp = <D(ap)/p di3p sm i/^p, 



(9) 



while t/p = 0. The scattering angle ij/^ is related to the phase angle 
by Qfp = n - ij/p. The total Stokes Q and U parameters measured 
by the distant observer in the observer's frame is obtained by 
transforming the reference frame and integrating over the visible 
surface of the star (as seen from the planet): 



visible 



(10) 



The Mueller rotation matrix L is defined by (e.g. lStenfloll 19941) 



cos 2ySp 
- sin 2y6p 



sin 2y6p 
cos 2/3p 



(11) 



where /3p is the rotation angle to transform the scattering plane 
reference system into the observer's frame. It depends on the 
stellar surface element. Note that Stokes V is zero in both refer- 
ence frames because we assumed that the light emitted from the 
star is unpolarized. The polarization induced by the unresolved 
extrasolar planet is then determined hy Q/F and U/F. We note 
that non-Rayleigh t ypes of sc attering influence the shape of the 
polarization curves (ISeager e t al. 2000). We will consider these 
cases in a forthcoming paper. 

2.3. Rayleigh scattering 

Let us now consider the case when both the intensity and po- 
larization of the light scattered in the planetary atmosphere are 
distributed according to the Rayleigh law. This implies solving 
a self-consistent radiative transfer problem. Here we solved this 
problem under the assumptions (in addition to those mentioned 
in Sect.lTTT) that 

- the atmosphere is plane-parallel and static and consists of 
homogeneous layers, 

- the planet is spherically symmetric, 

- the radiation entering the atmosphere is unidirectional and 
unpolarized, 

- an incoming photon is either absorbed or scattered in the 
continuum, 

- an absorbed photon does not alter the atmosphere (model 
atmosphere includes effects of iiTadiation), 

- single scattering approximation (a photon might still be ex- 
tinct after scattering). 
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Thus we solve a standard radiative transfer equation for scattered 
polarized radiation of a given frequency 



dl 

dr 

with the total source function 



S = 



KcB + (TrSr 



(12) 



(13) 



where and ctr are opacities due to continuum absorption and 
Rayleigh scattering, while 5r and B are the Rayleigh source 
function and the unpolarized thermal emission, respectively. 

The Rayleigh source function 5r is expressed as usual via 
the Rayleigh phase matrix P(fi, jj'), depending on the directions 
of the incident and scattered light. It has contributions from scat- 
tering of both the intrinsic thermal emission and incident stellar 
light. However, in the blue part of the spectrum, where the po- 
larization signal due to Rayleigh scattering is the largest, the in- 
tensity of the thermal emission is negligible compared to that 
of the reflected light. This is in contrast to the assumption by 
ISenguptal (120081) who neglected the polarization of the reflected 
light and took only into account the polarization of the thermal 
emission to interpret measurements in the B passband. In this 
case he had to also assume a non-negligible oblateness of the 
planet. This approach can only be justified to explain polariza- 
tion in red/infrared wavelengths which has not yet been detected. 

Taking into account that the incident stellar radiation is uni- 
directional (in case of a finite-size star the radiation comes from 
a given stellar pixel), we obtain 



Sj,(T,fl')^^P(M,fl')I(T,ti)^ 
4-7T 

where for single scattering 
/(r,ju) = /*e"^ , 



(14) 



(15) 



and T is the optical depth of the layer in the planetary atmosphere 
from T = at the top. In the same way the scattered radiation is 
reduced due to the optical depth on the way out of the atmo- 
sphere. 

The Stokes vector of the reflected light is obtained by solv- 
ing the radiative transfer problem as described above for a given 
vertical distribution of the temperature and opacity in a planetary 
atmosphere. The radiation flux is then obtained by integrating 
the Stokes vector over the illuminated planetary surface with a 
coordinate grid on the planetary surface. Our tests have shown 
that a grid of 6°x6° is sufficient to achieve the necessary accu- 
racy. Thus the phase function O(q') is in this case numerically 
evaluated. 

Our model includes the following opacity sources: (i) 
Rayleigh scattering on H, H2, He, H2O, CH4, and Thomson scat- 
tering on electrons, with all scattering species contributing to 
polarization, and (ii) continuum absorption due to H, H , H2^, 
H2 , He, He , Si, Mg, and Fe. The number densities of the 
species w ere calculated with a ch emical equilibrium code de- 
scribed in lBerd vugina et all ( 1200 3'). For an example calculation 
pr esented in Se ct.|24jwe employed standard model atmospheres 
by lAUard et al.l (1200 1ft . For interpreting real data, models of ir- 
radiated atmospheres are to be used. 

2.4. A comparison between the two cases 

The difference between the above two cases is in the angu- 
lar distribution of the intensity of the scattered light: in the 
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Fig. 1. Normalized phase functions for the Lambert sphere 
(dashed line) and the Rayleigh law (solid line). The calculation 
is for a circular planetary orbit with the inclination / = 30°, and 
in case of Rayleigh scattering for the wavelength of 4500 A. The 
differences occuiTing in the gradients of the two functions are 
caused by different angular dependences of the Rayleigh scatter- 
ing law and a Lambert sphere. For higher inclinations, the dif- 
ference in gradients is slightly modified. The dependence on the 
wavelength is a second order effect. The two curves have been 
renormalized to one at the maximum to make the comparison of 
the shapes more apparent. 



first case, it is assumed to be according to the Lambert sphere 
law, while in the second case, it is according to the Rayleigh 
law. Consequently the corresponding phase functions ^l(Q') and 
4>r(q') are certainly different, and it is important to evaluate the 
resulting difference in polarization curves. 

The shapes of the functions Ol(q') and Or (a) (normalized 
to their maxima) are shown in Fig. [T] Their behavior with the 
orbital phase (which differs from the phase angle a) is as ex- 
pected: the maximum illuminated area is observed at the exte- 
rior conjunction ("full moon", phase 0.5), while the minimum is 
at the interior one ("new moon", phase 0.0). Some differences 
between the two functions are seen in their gradients, which are 
due to the different angular dependencies of Rayleigh scattering 
and a Lambert sphere. The Rayleigh phase function also depends 
on the wavelength, but this effect is significantly smaller. Note 
that we normalized the plotted phase functions to one at their 
maxima to make the comparison of the shapes more apparent. 
In reality the two curves differ somewhat in scale depending on 
the detailed properties of the model atmosphere employed in the 
Rayleigh scattering case, with negligible influence on polariza- 
tion. 

The polarization curves resulting from these phase functions 
are shown in Fig.|2] The overall scale of polarization obviously 
differs for the two cases. Therefore we have changed the scale 
of the Rayleigh law results for an easy comparison of the curve 
shapes. The shapes of the polarization curves do not differ by 
a measurable amount for the case considered in this paper, i.e., 
when the reflected light is significantly dimmer than the direct 
stellar light and the planet remains spatially unresolved. This is 
because of the normalization of the Stokes parameters to the to- 
tal flux from the system. For resolved planets, the differences 
between the shapes of the polarization curves of the Rayleigh 
and Lambert cases increase but still remain small. This renders 
the Lambert sphere law a good approximation for the intensity 
of the reflected light, especially when obtaining orbital parame- 
ters of planets, which is the primary goal of our case study. To 
determine the composition and thermodynamic properties of the 
planetary atmosphere, a self-consistent polarized radiative trans- 
fer calculation as presented in Sect. l2.3l is required. 
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Fig, 2, Polarization curves obtained with the phase functions pre- 
sented in Fig. [T] The scale of the Rayleigh law curves has been 
adjusted to fit the Lambert sphere curves. The differences be- 
tween the shapes of Stokes Q and U are negligible given the 
current best accuracy of polarimetric measurements. Thus the 
Lambert sphere approximation for the intensity of the reflected 
light is justified for this case. 



2.5. Host star as a point-source 

We will compare the model described above for an extended host 
star with the simplified case, in which the star is treated as a 
point source. This case is applicable if the semi-major axis of 
the planetary orbit exceeds the radius of the star, so that all light 
rays from the star can be assumed to aiTive parallel on the planet. 
Hence the total flux at the observer with the normaUzation con- 
dition in Eq. (O is given by 



1 7 

1 pr 



where D is the distance between the star and the planet. The two 
parts in Eq. ( fT6b correspond to the direct stellar flux and the flux 
scattered by the planet, respectively. Since the incident radia- 
tion arrives parallel on the planet, it is sufficient to consider only 
one phase angle a. By comparing Eq. (fTSI l with the extended 
star model we find that D does not appear explicitly in Eq. ([8]l 
because it is implicitly contained there within the solid angle 
d£?p and in general differs for every individual surface element. 
In the limit of very large planetary orbits, i.e. approaching the 
point-source model, the distances between the planet and differ- 
ent stellar surface elements approach the same value, namely D, 
and can be explicitly expressed in Eq. ( fT6b . 

With respect to the scattering plane Stokes U' -Q and 



Q' 



pr 

cPD^ 



<l)(a) sin , 



which can be converted to the observer's frame with the Mueller 
rotation matrix L given in Eq. (fTTI) . 

3. Parameter study of polarization phase curves 

We perform in this section a parameter study and discuss the 
influence of different parameters on the QjF and UIF phase 
curves and on the degree of polarization P - ^JQ^ + U^j F. Since 
we showed in Sect. l2.4l that polarization curves calculated with 
the Lambert-Rayleigh approximation do not differ from those 
based on the self-consistent radiative transfer case for unresolved 
systems, here we employ the first case for simplicity. 

We considered the inclination of the planetary orbit /, the 
position angle of the ascending node Q, the eccentricity e, and 
the longitude of the periastron co. In addition we compared mod- 
els with extended and point-like stars to determine the parameter 
range for which the latter model serves as a good approximation. 

Most features of the polarized phase curves can be under- 
stood in terms of only a few quantities (cf. Eq. (fTTI)): 

L The polarization depends on the scattering angle ij/: it is the 
largest for 90° scattering at elongation and drops to zero for 
backward and forward scattering due to geometrical reasons. 

2. The greater the phase function O, the larger the polarization, 
since O describes the fraction of the incident light scattered 
towards the observer. The phase function reaches the maxi- 
mum at the smallest phase angle when the planet is farthest 
away from the observer and the minimum at the largest phase 
angle when the planet is closest to the observer. 

3. The polarization scales with the inverse square of the dis- 
tance D from the star to the planet, in accordance with the 
incident flux. In the case of elliptical orbits D varies with the 
orbital phase. 

4. The larger the radius of the planet, the larger the scattering 
surface and polarization. 

5. The linear polarization is always directed perpendicular to 
the scattering plane, i.e. perpendicular to the direction from 
the star to the planet as seen in projection on the sky (we 
neglect here effects of magnetic fields). Therefore, the linear 
polarization generally rotates while the planet orbits the star, 
so that Stokes Q and U are exchanged and change their signs. 

Both the radius of the planet and the semi-major axis of the 
planetary orbit only scale the amplitude of the polarized phase 
curves, and we will not further discuss them in this section. 
Uncertainties in the planetary radius do not influence the orbital 

(16) solution. We assume a - 0.05 AU, r-\R} (Jupiter radius), and 
Ri,-\ Rq. All other points in the list, including also the orbital 
parameters Q., e, co, modify the shape of the fight curves. 

3.1. Inclination of the planetary orbit 

First we discuss the inclination / of the planetary orbit. It is de- 
fined in a way that the values between 0° and 90° coiTespond to 
the counter-clockwise rotation of the planet as seen in projection 
by the observer, while angles between 90° and 180° indicate the 
clockwise rotation. The orbit is seen edge-on at an inclination 
i=90°. 

Figure[3]illustrates how the inclination / affects the linear po- 
larization at different times of the planet period, for simplicity in 
the case of a circular orbit. Basically it dictates how strongly the 
extrema of the polarization curves are pronounced. As always 
for circular orbits, we defined the time zero to by the passage of 

(17) the planet through the largest phase angle, i.e when it is closest 
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Fig. 3. Polarization curve dependence on inclination of the plan- 
etary orbit. The inclinations cover / = 0° (solid), / = 30° (dotted), 
! = 60° (dashed), / = 90° (dash-dotted), and / = 180° (dash-triple- 
dotted). The curves for / = 0° and /= 180° coincide in the top and 
bottom panels and are shown with solid lines. Further parame- 
ters: a = 0.05 AU, Q = 90°, e = 0, r = 1 and 7J* = 1 Rq. The 
orbital phase zero is defined by the largest phase angle (and by a 
position angle of 180° in the special cases of / = 0° and i- 180°). 
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Fig. 4. Polarization curve dependence on the position angle of 
the ascending node Q. Different curves correspond to Q = 0° 
(solid), O = 10° (dotted), Q = 20° (dashed), O = 45° (dash- 
dotted), and Q = 90° (dash-triple-dotted). Further parameters: 
semi-major axis a - 0.05 AU, inclination / = 75°, eccentricity 
e-Q, radius of the planet r-l Rj, and a stellar radius R*-! Rq. 
The orbital phase zero is defined by the largest phase angle. 



to the observer. The orbital phase represents the fractional time 
passed within one orbit relative to fo. 

When looking face-on onto the orbit, i.e. with / = 0° or 
i = 180°, the degree of polarization remains constant because 
the planet is always seen under the same phase angle. Still the 
rotation of the linear polarization leads to a sine-like variation of 
Stokes Q and U with the same amplitudes, but a phase shift of 
90°. 

The closer to edge-on we observe the orbit, the closer Stokes 
Q, Stokes U, and polarization approach zero for the orbital phase 
0.0. This minimum in polarization is very broad because the 
phase function and sin if/ (forward scattering) are simultane- 
ously small. Interestingly, both Q/F and U/F aie usually very 
small near to, except for inclinations below 30°, a helpful fea- 
ture for the determination of the zero point of the polarization 
scales (cf. Sect. |4]i. The polarization maxima occur near maxi- 
mum elongations, where sin^ = 1^ but slightly shifted towards 
smaller phase angles (i.e. closer to orbital phase tp = 0.5), be- 
cause for those phases the illuminated fraction of the planetary 
disk increases for the observer towards (p = Q.5. 

For the specific choice of Q = 90° in Fig. |3]all Q/F curves are 
symmetric with respect to phase if = 0.5, while all U/F curves 
are antisymmetric. This indicates that the shape of 2/^^ is inde- 
pendent of the direction into which the planet moves within the 
orbit and thus identical for orbits with inclinations / and 90°-/. 
On the contrary, U/F changes its sign for opposite rotation di- 
rections (cf. the solid and dash-triple-dotted curves in Fig. [3]). 
Note however note that these features depend on O. 



We also found that the extrema shift with inclination, not 
so much in the total degree of polarization, but clearly in Q/F 
and U/F. If we look for example at the ; = 60° case (dashed 
curve. Fig. |3]l, we find that the two main extrema of U/F occur 
relatively close to phase (p-0.5, even closer than the two maxima 
of Q/F, which is obvious from geometrical considerations. As 
we will clarify in the subsection below these properties of Q/F 
and U/F can be exchanged by variations of the position angle of 
the ascending node D.. 

Our model will allow us to significantly improve the mass 
determination of extrasolar planets, for most of which only a 
lower limit M sin / is known. According to the standard model 
of planet formation we expect that the angle between the incli- 
nation of the planetary orbit and the spin axis of the host star is 
small as in the case of the solar system. Indeed, it was found to 
be less than 1 1° in the sample of three transiting exoplanets for 
which the alignment was determined (IWinn etal.ll2005l 120061: 
Wolf et al. 2007). Therefore the measurement of the polarization 
curves, which are very sensitive to the inclination of the orbit as 
discussed above, represents a powerful method to constrain also 
the inclination of the stellar spin axis to within a few degrees. 

3.2. Position angle of fhe ascending node 

A variation of the position angle of the ascending node Q. rotates 
the orbit around the line-of-sight and is analogous to the rotation 
of the reference frame for polarization. The degree of polariza- 
tion remains unaffected accordingly, but the plane of linear po- 
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Fig. 5. Polarization curve dependence on the eccentricity. 
Different curves correspond to e = 0.0 (solid), e - 0.2 (dotted), 
e - 0.4 (dashed), and e - 0.7 (dash-dotted). Further parameters: 
a = 0.05 AU, / = 75°, Q = 90°, w = 90°, r^lR,, and = 1 Rq. 
The orbital phase zero is defined by the largest phase angle, i.e. 
when the planet is closest to the observer, which in this particu- 
lar case coincides with the periastron epoch due to our particular 
choice of w = 90°. 
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Fig. 6. Polarization curve dependence on the periastron longi- 
tude for eccentric orbits. Different curves correspond to w = 0° 
(solid), w = 45° (dotted), w= 135° (dashed), and w = 225° (dash- 
dotted). Further parameters: a = 0.05 AU, / = 75°, Q = 90°, 
e - 0.4, r - I Rj, and = 1 Rq. The orbital phase zero is de- 
fined by the periastron passage. 



larization rotates, as shown in Fig. |4] While increasing Q from 
0° (solid) to 45° (dash-dotted) the shape of Q/F changes from 
being symmetric to antisymmetric (with respect to (f = 0.5) and 
effectively takes over the properties of the U/F curve. The oppo- 
site is true for U/ F. At Q = 90° (dash-triple-dotted) both curves 
again reach the same shape as for Q - 0°, but with opposite sign. 

The polarization curves contain one ambiguity, since they 
remain unaltered when shifting f2 by 180°. This is of course an 
inherent property of polarization, because Stokes Q and U are 
not subject to change under a rotation of the plane of polarization 
by 180°. 

3.3. Elliptic orbits 

In the case of eccentric orbits we need to distinguish two addi- 
tional parameters that affect the polarization curves in a different 
manner, namely the eccentricity e and the periastron longitude 

CJ. 

For eccentric orbits we choose the periastron passage as zero 
to, because the polarization peaks shift towards the periastron, 
and the radial velocity method well constrains the periastron 
epoch for elliptic orbits. To allow a comparison with circular 
orbits we selected u) = 90° in Fig.|5]so that the periastron epoch 
coincides with fo used in Sects. [ni and D. 21 where it was defined 
by the largest phase angle. 

The eccentricity e influences the polarization curves in two 
ways when the semi-major axis is kept fixed: it shifts and scales 
the extrema of the light curves (Fig. |5]l. The shift of the ex- 



trema results from the varying orbital velocity of the planet. With 
increasing eccentricity the extrema group closer to the perias- 
tron passage, although maximum polarization is generally not 
reached exactly at the periastron. At the same time most extrema 
become stronger, in the case shown in Fig. |5]even all extrema, 
because the planet approaches the star much closer at periastron 
for eccentric orbits (with fixed semi-major axis). Recall that the 
incident flux on the planet and thus the scattering polarization 
scale with the inverse square of the distance D between the star 
and the planet. 

Initially it might sound surprising that all extrema strengthen 
for elliptic orbits in Fig.|5] since for circular orbits the maximum 
polarization always occurs at phase angles smaller than 90°, 
which occur when the planet moves through the part of the orbit 
around the apastron for our specific choice of (w = 90°. However, 
for very eccentric orbits the phase angle decreases quickly after 
the periastron passage, and the phase function exceeds the value 
of 0.5 already when the planet moves through the descending 
node. This happens when D is still much smaller than the semi- 
major axis and long before the maximum elongation. Moreover 
the phase angles at which the extrema occur move closer to the 
periastron due to the strong dependence of polarization on the 
distance D. The distance dependence of the polarization also ex- 
plains why the extrema in Fig. |5] are strongly scaled up near the 
periastron passage. By comparing the circular orbit (solid) with 
the most eccentric case (dash-dotted) we find indeed that the in- 
conspicuous minimum of U/F at if = 0.2 in the circular case is 
shifted to ifi - 0.03 and scaled up by about a factor of 10, the 
maximum in Q/F at i^ = 0.3 (circular case) is increased by a fac- 
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Fig. 7. Influence of secondary eclipses. The curves correspond 
to three different inclinations: / = 40° (solid), which is too small 
for the occurrence of eclipses, / = 60° (dotted), and / = 80° 
(dashed). The additional parameters correspond to a hypothet- 
ical hot Jupiter orbiting a giant with a = 0.109 AU, Q = 100°, 
e - 0, r - 9 Rj, and a stellar radius = 13.3 Rq. The orbital 
phase zero is defined by the largest phase angle. For other stellar 
or planetary radii the curves are scaled in the degree of polariza- 
tion, but their shapes remain the same. 



tor of 5, and the maximum of U/F (p - OA (circular case) is 
scaled up only by about a factor 2. 

In Fig.|6]we demonstrate the effect of the longitude of the 
periastron w, which in principle rotates and positions the elliptic 
orbit within the orbital plane. This results in an enhancement of 
the polarization curve extrema that lie closest to the periastron 
epoch. Therefore some Stokes U curves in Fig. |6] exhibit pro- 
nounced negative peaks that are not apparent for other values of 
oj. On the other hand, a variation of co can not produce negative 
(or positive) polarization if it is absent in the corresponding cir- 
cular case with the same orbital orientation. Such an example is 
visible in Stokes Q of Fig.|6l where strong peaks with negative 
Q are nonexistent for any oj, as expected from the corresponding 
circular orbit (solid line, top panel, Fig.jSj- 

From Fig. |6] we find some symmetry properties that occur 
with a change of cj. The curves with oj = 90° +y (or co = 270° + 
with y any angle, are just mirrored at(p- 0.5 in the total degree 
of polarization and in Stokes Q, while in Stokes U they are mir- 
rored at ^ = 0.5 combined with a sign change (see dotted and 
dashed curves). In Stokes Q and U this symmetry relation de- 
pends on the particular choice of Q, but it holds for any Q. in the 
degree of polarization. 



3.4. Size of the host star and secondary eclipses 

Previous simulations of the polarization signal from exoplanets 
always assumed that the light from the parent star falls parallel 
onto the planet (even in the case of hot Jupiters), a model we will 
refer to as "point-like" in the following. However, the assump- 
tion of parallel rays incident on the planet clearly fails in the case 
of hot Jupiters orbiting giants. This necessitates a correct mod- 
eling with an extended, finite-sized host star, which takes into 
account the strong distance dependence of the resulting polar- 
ization. 

A first immediate effect of models with extended stars is the 
possible appearance of transits and secondary eclipses. During a 
secondary eclipse the planet moves behind the star as seen from 
the line-of-sight, so that Stokes Q and U drop to zero, assuming 
that the direct stellar light is unpolarized. This effect is illustrated 
in Fig.Qwhere the inclination is varied. A secondary eclipse oc- 
curs if the inclination lies in the range 90°-/e < / < 90°+/e, 
where 4 depends on the relative size of the star and the plane- 
tary orbit. The effect of the eclipse on the polarization curves can 
be significant and provides a very sensitive tool to diagnose the 
orbital inclination and the orbital phase. Note that the planetary 
transits lead to much smaller relative signatures in the intensity 
light curves because only the direct stellar flux is slightly re- 
duced. In Fig.|7]the transits remain invisible although they were 
modeled in our calculations. 

The direct comparison of the point-like and extended star 
models reveals additional differences of the curves apart from 
eclipses (Fig. [8]l. First we consider the maximum degree of po- 
larization Pmax reached during one orbit as a general measure to 
compare the two models (Fig. |8^). The most significant devia- 
tions naturally appear for small planetary orbits, while the two 
models converge for large orbits as expected. We find that the 
two curves representing the two models cross at about a semi- 
major axis a/Ri, = 1.8. The exact location of this cross-over 
depends on the limb-darkening of the star, but always occurs 
around a/R-t. - 2. With a solar limb-darkening the cross-over 
happens at a//?* =2.3. 

For planetary orbits with semi-major orbits larger than 
a/R* ~ 2 the point-like star model always somewhat underes- 
timates the maximum polarization compared to the realistic ex- 
tended star model. The difference slightly depends on the stel- 
lar limb-darkening, but amounts to about 10% at the distance 
fl//?* ~ 3 and then continuously reduces for larger orbits to about 
5% at the distance of known hot Jupiters (a/Ri, ^ 10) and be- 
low 1% at a/R-^ ^ 60 (almost 0.3 AU for solar analogs). These 
deviations arise because the solid angle occupied by the star (as 
seen from the planet) is properly calculated in the extended star 
model, which leads to a slightly increased incident flux at large 
distances. 

For very small planetary orbits (a//?* <2) the point-like star 
model leads to a maximum polarization that is much too high, up 
to a factor of two higher than the realistic extended star model. 
At a first glance it might sound surprising that the extended 
star model results in a smaller maximum degree of polarization 
despite the planet being much closer to the stellar surface. At 
a/Ri, -2 the distance planet-stellar surface in the extended star 
model is halved compared to the point-like model, which in prin- 
ciples multiplies the incident flux by a factor of four. However, 
at these short distances only a part of the stellar hemisphere is 
visible from the planet, and the effective distance to points close 
to the limb is still relatively large, which almost compensates the 
proximity to the disk center, although the incident flux stays a bit 
stronger in the extended star model. 
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Fig. 8. Comparison of models with extended (solid) and point-like (dotted) stars, a) The panel illustrates the maximum degree 
polarization during one planetary orbit as a function of the semi-major axis a, which is given in units of the stellar radii. The relative 
difference is given in the panel underneath, b) Stokes Q curve at a distance a//?* = 1.4 and relative difference between the extended 
(solid) and point-like (dotted) stellar models. Note that Stokes U remains zero for the chosen parameters, c) Same as panel b but at 
a distance a/R-i, = 3.0. The following parameters were chosen for the calculations: / = 90°, f2 = 90°, e = 0, r = 1 and, in the models 
with a finite sized star (solid), ^* = 1 ^o- The dashed lines indicate the zero level. 



The crucial point is that in the extended star model the in- 
cident radiation does not arrive anymore parallel to the planet, 
which results in contributions to both Stokes Q and U and par- 
tial cancellation of polarization, while in the point-like model 
all incident (parallel) rays add up to the same direction of linear 
polarization. Assuming that light from the disk center gives rise 
to positive Stokes Q and zero U for the distant observer, then 
light from the limb starts to contribute mainly to Stokes U if the 
scattering plane is rotated by more than 22.5°, although the to- 
tal Stokes U cancels out due to symmetry reasons. If the planet 
orbits very closely to the stellar surface, some incident radiation 
would even arrive with angles greater than 45° relative to the 
stellar disk center, which would in turn lead to contributions to 
negative Q. At a//?* = 2 the star has an angular size of about 25° 
viewed from the planet. This is just the limit above which con- 
tributions to positive Stokes Q start to decline significantly. It is 
thus not surprising that the degree of polarization in the extended 
star model is actually smaller if the planet moves very closely to 
the stellar surface. 

Let us now look at the differences of the shape of the po- 
larization curves resulting from the point-like and extended star 
models, respectively (Fig. [St&c). For large orbits, i.e. beyond 
the cross-over distance of the two models, the maxima of the 
two light curves occur at about the same phase angles (Fig. [8};). 
The comparison of the maximum degree of polarization Pmax 
shown in Fig. [8^ is thus sufficient to estimate the deviations of 
the two models: the polarization curve is just scaled according 
to the relative difference of P,mx, except that we have to con- 
sider the possible secondary eclipses in the extended model. For 
very small orbits, i.e. inside the cross-over distance, the shape 
of the light curves start to differ more significantly between the 



two models (Fig. [S}?). The maxima of the light curve start to 
shift toward maximum elongation in the extended model (solid 
curve). In general the position of the polarization maxima is de- 
termined by optimizing the effect of two factors, the increasing 
phase function (Eq. Q) and the optimal scattering angle lead- 
ing to particularly large polarization, which is 90° for Rayleigh 
scattering. In the extended star model with very small orbits the 
planet is illuminated from different directions, giving contribu- 
tions to the scattered light with various phase and scattering an- 
gles. Contributions from stellar surface elements that lie close to 
the observer gain particular importance because the correspond- 
ing phase function is largest compared to the rest of the stellar 
surface. This effectively shifts the polarization maxima towards 
smaller phase angles. The polarization curves can easily differ 
by 20^0% (Fig. [8J3). In fact we find that the comparison of 
^'max (Fig. [Hh) underestimates the relative errors introduced by 
the point-like star model if a//?* <2. 

4. Fitting procedure to simulated data 

To derive the orbital parameters from the polarization data of 
extrasolar planets we performed a ;^f^-minimization. With the 
multitude of different Stokes Q and U phase curves presented 
in Sect. [3] one can wonder to what extent they are unique and 
how easily the correct parameters can be retrieved. In this sec- 
tion we test the robustness of the fitting procedure. We consider 
two cases of circular and elliptical orbits with simulated data and 
try to recover the original parameters with;t'^-minimization. 

The orbital phases of the simulated data were randomly cho- 
sen within 10 bins of the orbital period to mimic the situation 
with real observations in which one would try to obtain a full 
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phase coverage without actually achieving identical time steps. 
For these phases and a set of selected orbital and planetary pa- 
rameters we computed Stokes QIF and UjF, which defined our 
input or true polarization curves. Then we applied a random 
Gaussian noise to the true curves resulting in a set of simulated 
data with different signal-to-noise levels in the two studied cases. 
For a given signal-to-noise ratio we defined the signal to be the 
difference between the maximum and minimum of the curve ei- 
ther in Stokes Q or U, whichever was greater, while the noise 
level corresponded to one standard deviation of the Gaussian dis- 
tribution. 

Error bars of the best-fit parameters were evaluated with 
the help of Monte Carlo simulations. In the case of simulated 
data we have the advantage to know the "true" solution which 
is given by the chosen input parameters. For both considered 
cases we calculated 200 samples of measurements by applying 
Gaussian noise to this "true" solution in the same way as de- 
scribed above for the simulated data. In fact, the simulated ob- 
servations (Figs. |9] and [10] panels a and b) can be considered 
as particular examples of such Monte Carlo samples. The x^- 
minimization was applied to each Monte Carlo sample, resulting 
in 200 sets of best-fit parameters. The standard deviations from 
the true values define the Icr error bars for every free parameter. 

In the first simulated case (Fig. |9]) we assumed an elliptic 
orbit, given by 100 data points, and applied a moderate signal- 
to-noise ratio S /N -5 to the original curves (solid, panels a and 
b). In the fitting procedure we simultaneously optimized seven 
parameters: the planetary radius r, the orbital inclination /, the 
position angle of the ascending node Q, the eccentricity e, the 
longitude of the periastron oj, zero point shifts of Stokes Q and 
U. These shifts in the polarizations scales can account for back- 
ground polarization either due to circumstellar dust or interstel- 
lar scattering, although the most ideal situation is of course to 
measure possible constant off-sets. 

The Stokes Q and U phase curves are very well reproduced 
by the fitting procedure (Fig.|9^&b, dashed lines). All seven pa- 
rameters returned by the ^^-minimization agree well with the 
original input parameters (Table [1). In panels c-h of Fig. |9] 
we show two-dimensional projections of in the parameter 
space. The minimization procedure resulted in a minimum value 
^min ~ ^ degree of freedom. There exists however a sec- 

ond identical minimum, not indicated in the figure, in which O 
is shifted by 180°. Apart from this ambiguity, which is discussed 
in Sect. 13.21 the minimum of x^ is unique and well constrained. 
We found that the moderate signal-to-noise S /N = 5 is already 
sufficient to determine the true parameters within a few percent. 

The best-fit solutions of the Monte Carlo samples (dots in 
panels c to h. Fig. |9]l are of course centered around the true so- 
lution, which explains the small relative shift of the;t'^-contours, 
which indicate the best-fit solutions for the simulated data. As 
statistically expected this shift is smaller than the Icr error. 

In the second test case (Fig. [TOl i we considered only 50, 
i.e. half as many, data points and we considerably reduced the 
signal-to-n oise ratio down to 1.5, wh ich is more representative of 
the data bv lBerdvueina et an(l2008h . Here we assumed a circular 
orbit so that the longitude of the periastron becomes undefined. 
For the ;t^^-minimization we nonetheless kept the eccentricity e 
as a free parameter, but we fixed w to 90°. For very small eccen- 
tricities o) is not well constrained by the data but has also negli- 
gible influence on the polarization curves. Therefore, we are left 
with six free parameters. 

Even under these more difficult circumstances the fitting pro- 
cedure proved to be very robust. The original curves are again 
quite well reproduced (Fig. [TOh&b) and the input parameters 



well identified (Table [TJ. The ^-minimum Ot'^m - 119 per de- 
gree of freedom) is still well defined but not as pronounced and 
broader than in the first test case, which increases the uncertain- 
ties in the fitted parameters. Note the ambiguity in O, which is 
visible in panels c and e (Fig.fTOli. We conclude that the data of 
the second test case are still good enough to accurately approxi- 
mate the orbital parameters despite the small signal-to-noise. 



5. Conclusions 

We developed and tested the tools for the interpretation of po- 
larimetric data from extrasolar planets. We employed (i) an an- 
alytical approach based on the Lambert sphere approximation 
for intensity combined with Rayleigh scattering for polariza- 
tion and (ii) a numerical method for self-consistent polarized 
radiative transfer in a scattering and absorbing atmosphere. We 
demonstrated that these two approaches lead to somewhat dif- 
ferent intensity phase functions but negligibly different shapes 
of polarization curves. We focused on the influence of various 
orbital parameters on the Stokes Q and U curves for the gen- 
eral case of elliptic planetary orbits and carried out a parameter 
study. The parameters include the inclination, longitude of the 
ascending node, eccentricity and the periastron longitude of the 
planetary orbit. The shapes of the polarization curves are not 
only very sensitive to changes in all of these parameters, but are 
also unique for every set of parameters except for a 180° ambi- 
guity in Q that results from the intrinsic properties of Stokes Q 
and U. 

We found that the finite size of the parent star should be cor- 
rectly taken into account for very small orbits. Otherwise the 
approximation of a point-like star would lead to errors in model- 
ing polarization exceeding 10% if the planet resides closer than 
about three stellar radii from the surface of the star. For known 
hot Jupiters this error is reduced to about 5%, which might be 
acceptable depending on the polarimetric precision of the data. 

Fits to data with a;^'^-minimization proved to be very stable 
even for data with a signal-to-noise ratio close to one if a suffi- 
cient number of data points (typically 50 or more) are available. 
The Stokes Q and U curves are sufficiently unique that the fitting 
procedure can reliably determine up to seven free parameters si- 
multaneously. 

Polarimetric studies of extrasolar planets open up new ex- 
citing possibilities to explore these remote worlds, even when 
considering merely orbital parameters. Our analysis can derive 
the inclination of the orbit for example, which improves the de- 
termination of the planetary masses, for which currently often 
only a lower limit M sin / is known. The Stokes Q and U curves 
also constrain Q, the position angle of the ascending node of 
the planetary orbit. The orbital parameters of the planet could 
also lead to better estimates of the orientation of the stellar spin 
axis, which appears to be well aligned with planetary orbits. In 
possible transiting systems of hot Jupiters orbiting giant stars 
secondary eclipses can lead to drops in polarization on the or- 
der of 50% of the total signal, compared to a reduction in the 
flux of about 1% during transits. The most significant advances 
in the study of extrasolar planets will result however from future 
observations with greater polarimetric sensitivity that will allow 
us to study directly the chemical and physical properties of the 
planetary atmospheres. 
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Table 1. Parameters of the models used to test the fitting procedure with simulated data. 



Case 


S/N 


a 
[AU] 


r 


n 


n 
n 


e 


0) 

n 


Shift Q/F 
[xlO-*] 


Shift U/F 
[xlO-*] 


1 (Fig.ig 


5.0 


0.03^ 
0.03~ 


1.5 
1.5 ±0.1 


60 
61 ±3 


120 
124 ±4 


0.2 
0.19 ±0.03 


200 
199 ±4 


-2.0 
-2.1 ±0.1 


-0.7 
-0.6 ±0.1 


2 (Fig.m 


1.5 


0.03° 
0.03° 


1.5 
1.4 ±0.2 


75 
80 ±9 


90 
94 ±8 


0.0 
0.0* 


c 

c 


1.0 

0.9 ±0.3 


-0.5 
-0.4 ± 0.3 



° Fixed parameter 

* The uncertainty of e is 0.04 
Undefined parameter; fixed to 90° in the fitting procedure for eccentric orbits 
In each case the top row refers to the true parameters employed to simulate the data, while the bottom row gives the result obtained from the 
^^-minimization. The Icr errors correspond to the standard deviation of each parameter as obtained from the Monte Carlo simulations. 
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Fig. 9. Results of the fitting procedure using simulated observations with S/N = 5. a) and b) Input curves (solid) for Stokes Q 
(top) and U (bottom), from which the simulated data (dots) were generated by adding Gaussian noise with S /N-5, and the curves 
resulting from the best fit (dashed) to these data. Input and best-fit parameters are listed in Table[T](case 1). The orbital phase zero is 
defined by the periastron passage, c) to h);^'^ -contours (solid) as functions of the seven free parameters illustrated by a representative 
selection of two-dimensional projections within the parameter space. The contours show the differences to the minimum value for 
the levels Ax^ = 1.00, 2.71, and 6.63, whose projections to the individual parameter axes correspond to the 68.3%, 90.0%, and 
99.0% confidence intervals, respectively. The Monte Carlo sample solutions are illustrated with dots. [See the electronic edition of 
the Journal for a color version of this figure.] 
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